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Abstract  Transport  theor>  based  on  the  relaxation  time  formalism  has  been  applied  to  bismuth;  the  results 
are  used  in  Part  II  of  this  paper  to  determine  the  figure  of  merit  of  Bi  in  energy  conversion  processes.  Using 
the  Jones  Shoenberg  model  for  bismuth,  analytic  expressions  have  been  derived  for  the  electrical  resistivity, 
thermal  conductivity,  and  for  the  Hall,  Seebeck  and  Nernst  coefficients.  The  Boltzmann  transport  equation  was 
solved  for  the  perturbed  distribution  function  using  anisotropic  relaxation  times.  The  result  was  then  introduced 
in  the  transport  integrals  for  the  electric  current  and  for  energy'  flux  to  obtain  the  phenomenological  equations 
for  each  set  of  charge  carriers  associated  with  a  given  ellipsoid.  The  contributions  of  each  group  of  carriers 
were  then  added  in  the  common  symmetry  coordinate  system  of  the  crystal  to  obtain  the  above-mentioned 
transport  coefficients.  To  derive  analytic  expressions,  it  was  necessary  to  consider  the  special  cases  where  the 
magnetic  field  is  aligned  with  each  of  the  three  symmetry  axes  and  to  pass  to  the  limit  of  very'  low  or  very  high 
magnetic  fields. 


i.  PRELIMINARIES 

Recently  considerable  interest  has  developed  in  the  application  of  the  Nernst  and 
Nernst- Ettingshausen  effects  to  energy  conversion;  this,  in  turn,  has  stimulated  a  search 
for  suitable  device  materials.  Among  the  earliest  experimental  studies  in  this  direction 
were  measurements  on  the  galvano-thermomagnetic  properties,  figure  of  merit,  and  general 
device  performance  of  Bi  and  Bi  Sb  alloys  [14].  These  experiments  provide  an  opportunity 
for  checking  out  numerical  predictions  based  on  transport  theory;  if  the  answers  are  in 
agreement  with  experiment,  then  the  theory  can  be  used  in  further  investigations  con¬ 
cerning  optimal  operating  conditions  of  the  device.  The  purpose  of  our  work  is  thus  two¬ 
fold.  In  Part  I  transport  theory  is  developed  in  some  detail  to  establish  several  new 
features;  namely:  (a)  obtaining  the  solution  of  the  Boltzmann  transport  equation  for 
the  very  general  case  of  anisotropic  relaxation  times,  arbitrary  magnetic  field  strengths, 
and  temperature  gradients,  (b)  deriving  general  thermodynamic  relations  for  the  over-all 
Seebeck  and  thermal  conductivity  tensors  in  terms  of  one-band  contributions,  (c)  intro¬ 
ducing  anisotropic  mobilities  into  the  equations  of  interest,  (d)  establishing  the  necessary 
refinements  of  earlier  transport  theories  [5,6]  which  bring  theory  in  accord  with  experi¬ 
ments,  and  (e)  setting  up  all  explicit  formulae  needed  for  later  use.  In  Part  II  we  check  the 
theoretical  predictions  against  available  experimental  data  for  Bi  and  then  utilize  the 
theory  for  calculating  the  appropriate  figures  of  merit. 

The  band  model  which  is  used  will  be  introduced  in  Section  3;  however,  it  is  well  to 
remark  here  that  vve  shall  generalize  the  Abeles-Meibcxjm  [5]  treatment  of  Bi  by  taking 
into  account  the  Jones-Shoenberg  band  model  refinement  [7].  Also,  from  general 
considerations  it  emerges  [8]  that  materials  best  suited  for  energy  conversion  processes 
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based  on  the  Nernst  and  Nernst-Ettingshausen  effects  arc  intrinsic;  we  therefore  specia¬ 
lize  to  this  case  quite  early.  Other  simplifications  arc  introduced  as  needed  to  obtain  the 
final  results  in  tractable  form. 

We  begin  with  the  solution  of  the  Boltzmann  equation  in  the  relaxation  time  formal¬ 
ism.  The  procedure  parallels  that  of  an  earlier  derivation  [9]  to  which  it  reduces  for  the 
isotropic  case.  The  reader  is  referred  to  this  discussion  for  certain  details  and  symbols. 

One  should  note  at  the  outset  that  the  external  force  acting  on  charge  carriers  in  a 
band  in  zero  magnetic  field  is  not  simply  the  applied  electric  field  E,  but,  more  generally, 
the  gradient  of  the  band  edge,  V T6B  [9].  With  this  modification,  the  Boltzmann  transport 
equation  for  isotropic  media  reads : 

(Ze/hm-JJZe)  +  (v  x  H)/c].Vk/+  v.  Vr/=  -(/-/<,) f1,  (1.1) 

where  Z  =  + 1  for  holes  and  - 1  for  electrons,  v  =  8  and  SB  are  the  energies 

and  band  edge  energies,/ or /0  are  the  actual  or  equilibrium  distribution  functions  and  x 
is  the  relaxation  time;  the  remaining  symbols  have  their  conventional  significance.  We 
obtain  an  approximate  solution  to  equation  (1.1)  by  the  usual  method  [10]  of  setting 
f  ~fo  for  first  and  third  terms  on  the  left  and  /  =  /0  —  v .  ^(cf0/de)  in  the  remainder. 
The  first  objective  in  the  general  derivation  is  to  find  an  expression  for  the  quantity  *F. 
For  this  purpose,  one  proceeds  essentially  as  in  [9].  The  only  modification  required  to 
adapt  the  result  to  anisotropic  media  consists  in  the  replacement  of  t~1  by  the  tensor 
x-1  =  w  The  generalization  of  equation  (1.8),  [9],  is  thus  given  as: 

P .  v  -  (Ze/hc)(\  xH).(T.  Vkv)  =  (x  ~ 1 .  *F) .  v,  (1.2) 

where 

P  s  ZeVr(SB/e)  -  TWr(jiB/T)  -  (e/T)VrT  (1.3a) 

*  ZeVr(C/e)  +  l(fis  -  s)/T]VrT.  (1.3b) 


In  the  above,  piB  =  £  —  SB,  and  s  =  6  —  <j  Bl  where  £  is  the  Fermi  level. 

In  the  principal  coordinate  system  of  a  particular  ellipsoid,  one  can  write 

Vkr/ff  =  Vk,Vk*/**, 


where  the  primed  coordinates  serve  as  a  reminder  that  this  special  coordinate  system  is 
being  used;  the  quantity  on  the  right  represents  a  diagonal  entry  of  the  reciprocal  mass 
tensor  cf  =  m~l.  On  inserting  this  result  in  equation  (1.2)  and  applying  the  triple  scalar 
product  rule  to  the  term  (v  x  H) .  PF  .*q),  one  obtains : 

{P  -  (Ze/c)[H  x  PF  .  q)]  -  (x  " 1 .  ¥)} .  v  =  0.  (1.4) 

As  has  been  discussed  in  connection  with  the  comparable  equation,  (1.11b)  of  [9],  equation 
(1.4)  can  only  be  satisfied  by  requiring  that  the  quantity  in  curly  braces  vanish.  In  addition, 
we  introduce  a  new  variable  H**,  defined  as  *F  s  Y.  'F*.  Equation  (1.4)  can  then  be  re¬ 
arranged  to  read : 


*F*  -  (Ze/c)[(x  .  *F*).1q]  x  H  =  P. 


(1.5) 
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To  solve  the  above,  we  operate  with  V.  on  the  left ;  this  is  to  be  followed  with  the  operation 
.  qf  from  the  right  and  finally,  by  the  operation  x  H  from  the  right.  One  then  obtains : 

[(t  X  H  -  (Ze/c)  {[?.  x  H}]  .tf)  x  H  =  [(Y.P).^]  x  H.  (1.6) 

The  next  step  consists  in  evaluating  the  quantity  sandwiched  between  Ze/c  and  xH 
in  the  second  term  of  equation  (1.6).  The  relaxation  time  formulation  of  the  Boltzmann 
transport  theory  rests  on  the  assumption  that  energies  are  conserved  in  the  collision 
between  charge  carriers,  and  that  the  velocity  is  randomized  during  such  encounters. 
Also,  consistent  with  the  above  assumption,  scattering  of  carriers  from  one  valley  to 
another  is  assumed  to  be  included  in  the  relaxation  time  t.  Herring  and  Vogt  in] 
have  shown  that  in  these  circumstances  Vis  a  diagonal  tensor  in  the  ellipsoidal  coordinate 
system.  Let  us  set : 


Lx' 
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i 

mx. 

0 
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T  = 
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v  o 

:  q  = 
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;  H  =  (Hx ,  Hy.,  H. ). 

(1.7) 
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4»*  =  ( 4/* ,  4/*.  T*) 

On  performing  the  indicated  operations  in  the  term  under  consideration,  one  arrives  at 
the  following  result : 

It  .  {[(t .  4'*) .  q]  x  H}]| ,  q  =  V*  x  (F.  H),  (1.8a) 

where 

CSSt:rV||»||||x-‘||,  (1.8b) 

in  which  || m  is  the  determinant  of  the  mass  tensor.  Equation  (1.8)  may  be  checked  by 
substituting  from  equation  (1.7)  on  both  sides,  thus  obtaining  an  identity.  We  now  sub¬ 
stitute  equation  (1.8)  into  equation  (1.6)  and  apply  the  triple  vector  product  rule  to  the 
resulting  middle  term,  (, Ze/c )[*F*  x  (C  H)]  x  H,  to  obtain: 

[(V.  M'*).  <f]  x  H  -  (Ze/c){{ H.  T*)(C  .  H)  -  [H  (C  .  H)]Y*}  =  [(V.  P).  cf]  x  H.  (1.9) 

Next,  return  to  equation  (1.5)  and  operate  on  both  sides  with  .  H(C  .  H)  from  the 
right.  By  the  triple  scalar  product  rule,  the  resulting  second  term  on  the  left  vanishes 
identically,  leaving: 

(4'* .  H)(C  .  H)  =  (P .  H)(C  ,  H).  (1.10) 

so  that  on  substituting  this  result  into  equation  (1.9),  we  finally  obtain : 

[(V.  *F*).(f]  x  H  =  (Ze/c)(P.U)(C.H)  -  (Ze/t)[H .  (C  .  H)]*F*  4-  [(t.P).q]  x  H.(l.ll) 

We  have  hereby  succeeded  in  reformulating  the  second  term  of  equation  (1.5)  in  a  manner 
that  allows  this  relation  to  be  solved  for  ¥*.  Inserting  equation  (1.11),  one  obtains: 

T  r  T*  =  1  •  [P  4-  fc/c)2(P.  H)(C  .  H)  +  (Ze/c)[(?.  P)  .?]  x  H} 

1  +  (e/c)2[ H  .(F.  H)] 


(1.12) 
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from  which  the  distribution  function / can  be  calculated.  If  the  relaxation  time  is  a  scalar,  if 
VrT  =  0,  and  for  homogeneous  materials  where  VT(SB/Ze)  =  E  when  VrT  =  0,  the  above 
reduces  to  the  relation  cited  by  Shibuya  [12]  and  extended  to  anisotropic  relaxation 
times  by  Bullis  [13]. 


2.  THE!  GENERAL  PHENOMENOLOGICAL  EQUATIONS 

The  second  task  consists  of  writing  down  the  phenomenological  equations,  based  on 
the  transport  of  electric  charge  and  of  “kinetic  energy”  £,  as  given  by: 

.1  =  (Ze/47t3)Jy/d3k  =  —  (Z<?/47r3)Jv(v .  '*')(<%/&)  d3k  (21a) 

Jff.  =  (l/4n3)JevT d3k  =  -(l/4n3)Mv.  'V)(df0/?E)  d3k.  (21b) 

In  the  above,  the  integral  involving  f0  vanishes  identically.  To  obtain  an  explicit  formula¬ 
tion,  we  substitute  for  4*  from  equation  (1.12),  utilizing  equations  (1.3b),  (1.7)  and  (1.8b). 
After  carrying  out  the  required  operations,  one  obtains  the  following  expression  for  the 
.V  component  of  *P : 

4VA'  =  xxP  x’  +  (Ze/c)2xx.xy.x:4Hx/wy.m:)(HxPx.  +  HyPy.  +  H-.P..)- 
~  ( Ze/c)(zxz .  Hy  P.  /m.-  -  zxzy.H.Py./my) 

=  [Zezx>  +  (Ze3/c2)(zx.TrzIJ1X'/m.my,)]Vx.(t;/e)-\- 

+  [(Ze3/c2)(zx  t,.  t .Hx  Hy./my.m..)  +  (e2/c)(zxzy  H.-/my-)]Vy  (£/e)+ 

+  [{Ze3/c2)(zx.zy.z:.H Z.H x./my.m.,)  -  (e2/c)(zx  z:  lIy!m.  )]V.  (C/e)  + 

+  [(/<b  -  m V[1  +  (e2/c2)(xy-z.  Hl'/my-m.)]VX’T+ 

+  [(/<b  -  e)/T](xx.zy./wy.)[(Ze/c)H:.  +  (e2/c2)(z:.Hx.Hy./m:.)]Vy.T+ 

+  [(/<b  -  c)/T](zx.z:./mM-Ze/c)Hy.  +  (e2/c2)(zy.HxJI;./my)]V:.T,  (22) 

where 

A'  =  1  4-  (e2/c2)[Ty>z:  Hl  /wY  nizr  +  zx>xAlylmx>m:>  +  tv*t y>Hj>/wx>wy>].  (2.3) 


Corresponding  relations  for  Wy-  and  HC  arc  obtained  by  cyclic  permutation  of  the 
component  subscript. 

Equation  (2.2)  and  its  y',  z*  analogues  must  now  be  substituted  in  equation  (2.1). 
In  this  connection,  it  becomes  convenient  to  introduce  the  following  transport  integrals 
(/',  v' = 


g<;V) 

L\k'y) 


1 

47T3 


V2,  X; 


A' 


CF 


d3k 


1 

47I3C 


1  t% 

my  A'  ci: 


d3k 


1 

4n3c2 


"bA'lxxzyz:,r!  1  Pf0 

7  - XL  U  K. 

wx-myw.-A  cc 


(24a) 


(2.4b) 


(2.4c) 
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It  is  expedient  to  reformulate  these  integrals  by  expressing  the  energy  e  of  charge 
carriers  relative  to  the  appropriate  band  edge  as : 

/:  =  (h2/2)(kl>/nix>  +  k2/wy>  +  kl/m 2.)  (2.5) 

in  the  principal  coordinate  system  of  an  ellipsoid.  We  can  then  introduce  a  new  variable, 
=  kl'/m}/,  with  which  (2.5)  may  be  rewritten  as: 

«  =  (/,72)ie  =  m2.  (2.6) 

A' 

The  above  equation  characterizes  a  set  of  concentric  spheres  in  ^-space.  Henceforth,  we 
assume  that  the  t;>  are  functions  of  r.  only.  In  this  event,  we  may  write  d3^  =  4tt£2  d£, 
whence : 


d3k  =  (i»v*i»r.iii;r01/2  d3^  =  (4n/hz)(2d)^2(m^mY4tK>)lf2  da 
By  spherical  symmetry,  c2.  =  (l/3)<;2,  then: 

(m>:/h2)(dE/dky)2  =  =  2u/3. 

We  now  use  equations  (2.7)  and  (2.8)  to  rewrite  (2.4)  as: 

2v/2K.;»(..m.-.)1/2  ”  T,y+1/2  3/0 

**  3jt2mi;./i3  „  A'  ft: 


w»v)  =  2v/2(h)v  i)))..»i-)1/2  r  vr/^%, 

3n2cm;/mv>h^  0  A'  ?/: 

2J2 


3n2c2(mxnny*m::»)ll2h3  J0 


A' 


ft 


(2.7) 

(2.8) 


(2.9a) 

(2.9b) 

(2.9c) 


Utilizing  either  equations  (2.4)  or  (2.9),  the  \*'  component  of  the  current  density  vector 
J*'  =  -(Ze/47r1)njr2.4/v(5/0/ft)d3k  becomes: 

J*'  =  Ze(Ze/C,1v''  +  ZeaLiHax.WAi/e)  +  Ze(e2G(?Y)Hz.  +  Ze'LJIJijV,  ,{C/e)+ 
+Zc(-t-2H,.Gf=')  +  Ze3LlHx.HeWAl M  +  (Ze/T)[(K(?%  -  K(p)  + 

+  c2H2AL1(,h  L,)]V,.T 


+  (e2/T)[//r.(G^>'V„  -  Gf  ‘">)  +  ZcHx.Hy.(Llfill  -  L2)]V,,T+ 

+  (e2/T)[-H>,(Gtr-')fin  -  G(2  :))  +  ZeHy  H:.(L ,/<„  -  L2)]V.T.  (2.10) 

The  corresponding  quantities  Jr  and  Jr  are  obtained  by  cyclic  permutation  of  the  co¬ 
ordinate  indices.  Furthermore,  one  can  show  that  Jq>  =  —  (1  /47T3) ^Iev1^  x\df0/de)  d3k  is 
specified  by  a  relation  like  equation  (2.10),  except  that  one  power  of  Ze  is  deleted  from 
the  right  and  that  the  subscripts  g,  k ,  and  /  in  equation  (2.9)  arc  raised  by  one  unit.  The 
remaining  components  are  again  found  by  permutation  of  coordinate  indices. 

The  information  discussed  above  can  be  succinctly  summarized  by  the  following 
relationship: 


J 

<T 

P 

6 

i 

W/e) 

_J  O'. 

u 

i 

i  Y  J 

VT 

(2.11a) 
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where  in  the  principal  coordinate  system  of  each  ellipsoid : 

~e2Kp  +  eiL1Hy 

Ze'Gp'H.. 

-Ze'Gp'H,? 

+  eiLlHxHy 

+  e*L,Hx.Ht. 

—ZeiG{*',~>Ht. 

e2K'p 

ZeWp'Hy 

+  c4L,K;(.//y. 

+  e*LlH2. 

+  e‘>LlHyH, 

Ze'Gp'Hy 

-Ze'G'p'Hy 

+  e2K{p 

.  +  e4L,  HXH, 

+  eiLlH./H,, 

+  e4Z.,//| 

~(Ze/T)(K[*>nB  -  *</  >)  + 
+(Ze3/T)(Ltp>„  -  L2)H2X- 

~(e2/T)(Gp\iB  -  G(2  ,  ))H.  + 
+(Ze3/T)(Llt i„  -  LJH..H,. 

+{e2/T)(G,*‘  'nB  -  Gf-'')H,  + 
_  +  (Ze3/T)(L,/is  -  L2)Hx.H: 


+(*J/T)(G,itVV*  -  &r  ')H,+ 

+  (Ze3/T)(L^a  ~  L2)Hx-H, 

(Ze/T){KYh‘a  ~  K'P)+ 

+{Ze3/T)(L^„  -  L2)H2 

-(e2/T)(G<r\tB  -  Gp  >)Hx.+ 
+  (Ze3/T)(L,n„  -  L2)HyH:. 


(-e2/T)(G'r  >HB  ~  Gf '•>)//,.+ 
+  (Ze3/T)(LlnB  -  L2)Hx  H,, 

(e2/T)(Gp  >VB  -  Gp  ')Hx.+ 

+  (Ze3/T)(LIftB  -  L2)H,.H, 

+  (Ze/T)(KpnB  -  K(p)  + 

+  (Ze3/T)(LlfiB  -  L2)H2 


(2.11c) 


ZeKp  +  Ze3L2H2x. 


U  = 


-e2Gp  'H, 

+  Ze3L2HxH , 


e2G{21)Hy. 

+  Ze3L2HxH. 


e2G2‘,)Ht, 

+  Ze3L2Hx.Hy 

ZeK'p 
+  Ze3L2H2 

-e2Gp,Hx 
+  Ze3  L2H  y  H . 


-e2G{p'Hy 
+  Ze3L2Hx.Ht. 


e2Gp>Hx. 

+  Ze3L2H,.Ht. 

ZeK'p 
+  Ze3L2H2 


(2.1  Id) 


and 

\i/T)(Kp%  -  Kp) 

+  (e2/T)(L2t iB  -  L3)H2x. 

(-Ze/T)(Gp‘>»B  -  Gp  >)H: 
+  (e2/T)(L2v„  -  L3)Hx.Hy. 

(Ze/T)(GpPa  -  Gp  ')H,. 
+  (e2/T)(L2t,B  -  L,)HX  H, 


(Ze/T)(Gp,)tiH  -  Gp  ')H 
+ie2/T)(L2n„  -  L3)Hx.Hr. 

(1  /T)(K'pt,B  -  K>p) 
+{e2/T)(L2vB  -  L3)H2 

-(Ze/T)(G'p  >t,B  -  G'p-  ')HX. 
+  (e2/T){L2»B  -  L3)Hy-H._. 


-(Ze/TWr  ^a  ~  G?V,)H, 

+  (e2/T){L2t,B  -  L%)Hx.Ht. 

(Ze/T)(G'p  '» B  -  Gp  ')Hx. 

+  (e2/T)(L2»„  -  L3)Ht  H;. 

( \/T)(K'pn„  -  K<p) 

+  (e2/T)(L2n„  -  L3)Hl 


(21  le) 


To  avoid  the  complexities  inherent  in  the  use  of  equation  (2.11),  we  now  make  another 
in  a  series  of  approximations  by  taking  each  t;.  to  be  independent  of  e.  In  this  event,  we 
can  rewrite  equation  (2.9a-c)  as : 

K™  =  (r,./m,.A')Sk 

G<'V)  =  (x;/Zy/mymv.A'c)Ss 

Li  =  (xxTy.T../mxmy.m:A'c2)S„ 


(2.12a) 

(2.12b) 

(2.12c) 
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in  which  a  new  transport  integral 


S„  = 


hil  K-»v»v)1/2  7  _» + 1/2 jfjn  j - 

3ti2  F  J  &  dC 


(2.13). 


has  been  introduced;  n  refers  to  the  index  k, g, l  in  equation  (2.12).  The  quantity  A'  is 
given  by  equation  (2.3). 

The  zero  magnetic  field  charge  carrier  conductivity  partial  mobility  components  are 


u>:  =  eT>:/mx. 


(2.14) 


in  the  principal  coordinate  system  of  each  ellipsoid.  This  relation  applies  since  we  assumed 
earlier  that  zr  does  not  depend  on  e;  moreover,  the  u,:  represent  mobilities  in  zero 
magnetic  field.  Equation  (2.3)  now  becomes: 


A'  =  1  +  uy’U2  Hl’/c2  +  uX'U2’Hy‘/c 2  -f  ux  uy  H2  /c 2, 


(2.15) 


and  the  matrices  (2.11)  reduce  to: 

a  =  S,eM 
S^ZQtM/T 
tf  =  ZS2  M 
V  =  Q2M/eT, 

where  Q{  =  StfiB  —  Sf+1,  and  where 


Ms  — 
A' 


/ 

[l  +  'w(~ )  J 

«w  ( 

Z-f  +  »,  [l  +  «x«v(~)  ] 

+  u,,(-z!k  +  Ux"-£.) 


(2.16a) 

(2.16b) 

(2.16c) 

(2.16d) 


/  Ht. 

u,u,(z!±  +  uW) 

4^)1 


.  (2.16e) 


The  above  relations  represent  the  end  result  of  the  statistical  transport  theory  for  a 
group  of  charge  carriers  associated  with  an  ellipsoid,  in  terms  of  the  principal  coordinate 
system  of  that  particular  ellipsoid.  We  now  wish  to  correlate  the  tensor  entries  with 
quantities  that  are  experimentally  determined.  For  this  purpose,  it  is  convenient  to 
introduce  phenomenological  equations  in  partially  inverted  form  [14],  namely: 


r  «-* 

P 

-i 

P 

j  ' 

.  T Js. 

ft 

—  K  . 

.VT. 

The  entries  in  the  above  relations  are  essentially  definitions.  Thus,  when  VT  =  0,  VC/e 
and  J  are  interconnected  by  the  resistivity  tensor,  p*;  when  J  =  0,  VC/e  and  VT  are  related 
by  the  Seebeck  tensor  Again,  for  J  =  0,  T3S  is  the  heat  flux  (energy  transfer  past  unit 
area  in  the  absence  of  a  net  particle  transfer),  which  is  related  to  the  temperature  gradient 
by  the  negative  thermal  conductivity  tensor,  —Ik.  ^Finally,  for  VT  =  0,  the  functional 
dependence  of  TJS  on  J  involves  the  Peltier  tensor,  11.  It  is  wrorth  noting  that  the  entries 
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in  the  latter  are  connected  to  the  |f  entries  by  the  Casimir -Onsager  reciprocity  condition 
niV(-H)  =  TpM( H);  furthermore,  />*/- H)  =  p}i{ H),  and  /clV(-H)  =  k^H). 

For  the  moment,  we  will  continue  to  consider  only  the  group  of  charge  carriers 
associated  with  a  given  ellipsoid.  In  this  event,  the  entries  in  equations  (2.17)  can  be 
expressed  in  terms  of  those  in  equation  (2.11a)  or  (2.16)  as  follows:  By  appropriate 
manipulation  of  the  results  in  [15],  it  can  be  shown  that: 


TJS  =  JQr  -  K/VT  -  (fiB/Ze) J.  (2.1 8a) 

Actually,  we  wish  to  ignore  the  contribution  of  the  lattice  thermal  conduction  processes 
to  and  define  a  new  quantity  TJy  by 


7  Jy  =  J. 

Applying  equation  (2.18b)  to  equation  (2.16),  we  can  then  write: 


4— > 

70] 

“  - 

.1 

S,<n\1 

Vl  M 
T 

v:  e 

7JV  +  ~r 

“M 

V7' 

Ze  J 

el 

(2.18b) 


(2.19a) 


(2.19b) 


Now  solve  for  VC/e  in  terms  of,]  and  VT  in  equation  (2.19a)  and  eliminate  VC^e  from 
equation  (2.19b):  also  replace  the  Q{  by  their  definition.  We  then  obtain: 


Comparison  between  equations  (2.20)  and  (2.17)  then  leads  to  the  following  identifications : 

^/>=  M-*<vy SfV,  (2.21a) 

where  it  is  clear  from  equation  (2.13)  that  Sx  depends  on  the  effective  mass  of  the  carriers 
under  consideration;  is  therefore  provided  with  the  index  v.  For  all  ellipsoids  lying 
within  a  given  band,  S(xv)  is  the  same:  this  is  also  true  of  the  mobilities  of  carriers  in  a  given 
band.  ^ 

One  can  determine  M“ 1(v)  in  equation  (2.21a)  from  equation  (2. 1 6e)  using  the  standard 


inversion  procedure.  This  leads  to  the  result: 

r  i  //<v> 

n(y)l 

• y  1 1  y 

A  v 

_z  - 

»vV  C 

C 

- .  V 

//<V'  1 

-7.W- 

1  p  S\"c 

/v  r  ' 

V 

c 

„  „  / 

7?’  1 

z>  Zy 

c 

c  u{e. 

(2.21b) 
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From  a  comparison  of  equations  (2.20)  and  (2.17),  the  general  Seebeck  tensor  is  found 
to  be: 


ft"  = 


ov 


(V) 


Z,eTS{? 


1  0  0 
0  1  0 
0  0  1 


ZjuT 


(») 


[1]  =  P("[1] 


(2.22) 


In  the  present  approximation,  therefore,  the  Secbcck  coefficient  for  a  single  valley  is 
isotropic. 

Finally,  the  electronic  contribution  to  the  thermal  conductivity  is: 


tMv)  _ 
K  />  — 


S(/'>S(,v)  _  s(f . 
~eTSfp 


Mw  = 


S(V)SM  _  S(v)-’ 

e2TS'/'2 


(2.23) 


Since  k^(v)  is  proportional  to^^  the  Wicdcmann-Franz  law  holds  in  the  approximation 
scheme  used  here. 

Having  disposed  of  the  one-group  problem,  it  remains  to  find  expressions  for  the 
totality  of  carriers  in  all  different  groups  of  relevance.  In  general,  the  physical  properties 
of  anisotropic  materials  are  characterized  in  terms  of  several  bands,  one  or  more  of  which 
is  of  the  multi-valley  form.  For  the  summation,  we  use  as  a  guiding  principle  the  fact 
that  the  individual  lluxcs  arc  additive,  whereas  the  gradients  VT  and  Vf/e  are  the  same 
for  each  group  of  carriers.  Furthermore,  one  should  recognize  that  equation  (2.17)  is  a 
thermodynamic  relation  that  takes  no  account  of  the  details  of  conduction  processes. 
Hence,  the  form  of  equation  (2.17)  remains  unaltered  for  the  totality  of  charge  carriers; 
J  and  Js,  then  refer  to  the  total  flux  vectors,  and  the  tensor  entries  specify  the  appropriate 
physical  properties  of  the  material  as  a  whole. 

Let  us  then  determine  the  total  current  density  from  equation  (2.19)  as: 

J(T)  =  X J(v)  =  X XV'VM<v).Vc/e  +  X (ZvQ(1v>M(v)/r)  VT 

V  V  V 

=  (£#"> )  •  VC/e  -  f  V  p(,Vv))  •  VT,  (2.24) 

x  V  /  V 


where  equation  (2.22)  was  used  to  arrive  at  the  second  line.  Likewise,  from  equations 
(2.20  2.23),  the  total  entropy  flux  (exclusive  of  the  contribution  of  the  lattice  to  the  heat 
conduction)  is  found  to  be: 

TJ(P  =  T  X  W  =  T  X  P(v)J(v)  -  X ^'(v)  •  (2.25) 

V  V  Y 

Now  eliminate  .J(v)  in  the  above,  using  equation  (2.19a):  then : 

T.][P  =  -  £k'<‘>  •  VT  +  X  Tp(,)(5*v)  •  VC/e  -  p(v,oJ<v)  •  VT).  (2.26) 

V  V 

For  VT=  0,  it  follows  from  equation  (2.24)  that  Jm  =  ■'VC/e  whence  the  total 
conductivity  and  resistivity  tensor  is:  v 


«  =  K  =  (l^v,r.  (2.27) 

According  to  equation  (2.24)  [see  also  equation  (2.17)],  when  J(r)  =  0,  the  total  Seebcek 
tensor  is  given  by : 

p  ss  V(C/e)/VT]j( T)=  0  =  p(v^v)). 


(2.28) 
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Finally,  the  total  thermal  conductivity  (exclusive  of  the  contribution  from  the  lattice)  is 
found  from  equation  (2.17)  as  —  T3(P/VT  when  J  —  0.  Using  this  latter  condition,  we 
can  eliminate  V£/c  between  equations  (2.24)  and  (2.26),  thereby  obtaining : 

+  T  V  (p(v))2<?(v)  -  t(2p(vWv)).P.(Zp(v)»w) 

V  V  '  V  '  '  V 

=  +  T  V  (p,v))2o(v|  -  t(Xp(,’S(v,)-K  (2.29) 

V  V  '  V 

Having  shown  how  to  determine  and'K'  from  one-group  contributions,  we  can 

now  apply  the  above  theory  to  the  case  of  Bi.  Due  to  severe  mathematical  complexities 
whieh  arise,  we  shall  not  explicitly  determine *k  according  to  equation  (2.29);  the  reader 
will  note,  however,  that  this  can  be  done  in  principle  through  a  knowledge  of  P(v)  and 
<r(v),  since  the  V(v)  occurring  in  equation  (2.29)  can  be  eliminated  via  the  Wiedcmann- 
Franz  formulation,  equation  (2.23). 


3.  TRANSPORT  TENSORS  FOR  THE  CONDUCTION  BAND  OF  BISMUTH: 

THE  TILTING  OPERATION 

Having  constructed  in  Section  2  tensors  for  various  transport  coefficients  in  the 
principal  coordinate  system  of  each  ellipsoid,  we  now  begin  the  task  of  deriving  the 
corresponding  quantities  for  Bi.  For  this  purpose,  it  is  necessary  to  consider  its  band 
structure. 

In  our  further  work,  we  shall  specify  the  orientation  of  ellipsoids  with  reference  to  a 
crystal  axis  system  in  which  the  binary,  bisectrix  and  trigonal  axes  constitute  the  x,  y, 
and  r  directions. 

Abeles  and  Meiboom  [5]  (AM)  assumed  a  conduction  band  exhibiting  six  whole  or 
six  half  ellipsoids  which  lie  along  the  threefold  binary  axes  and  hence  are  separated  by 
60°  in  the  x  v  plane.  They  also  postulated  a  valenee  band  in  which  two  half  or  two  whole 
ellipsoids  lie  along  the  r-axis.  The  final  formulae  are  independent  of  whether  the  bands 
are  assumed  to  contain  six  or  three  and  two  or  one  ellipsoids  respectively;  we  arbitrarily 
consider  the  ease  where  the  eonduetion  band  contains  three  ellipsoids  (labelled  v  =  1,2,3) 
and  the  valence  band,  one  ellipsoid  (labelled  v  =  4).  We  also  introduce  the  Jones-Shoen- 
berg  (JS)  refinement  [7],  in  which  the  conduction  band  ellipsoids  arc  rotated  about  their 
respective  binary  axes  in  such  a  manner  as  to  preserve  a  three-fold  symmetry  about  the 
trigonal  axis.  Also,  wc  extend  the  AM  treatment  to  the  Sccbcck  and  Nernst  tensors. 
In  recent  years,  experimental  evidence  for  two  additional  band  model  refinements  have 
appeared  in  the  literature;  namely,  (a)  the  non-parabolicity  of  the  bands  [16],  and  (b)  the 
occurrence  of  additional  bands  [17].  We  will  neglect  these  latter  complications  because 
introduction  of  (a)  in  the  theory  generally  leads  to  integrals  which  have  to  be  evaluated 
numerically,  and  both  (a)  and  (b)  lead  to  extremely  complicated  expressions. 

The  first  step  in  the  derivation  consists  in  introducing  a  similarity  transformation 
which  rotates  each  conduction  band  ellipsoid  about  the  appropriate  binary  axis,  (i.c., 
about  its  own  x'  axis)  by  an  angle  0(~  —6°),  thereby  bringing  its  x', y'  plane  into  proper 
alignment  with  the  x,  y  crystal  plane.  The  transformation  which  accomplishes  this  will  be 
denoted  as  the  “tilting  operation1’  and  the  conduction  band  ellipsoids  which  have  been 
subjected  to  this  operation  are  said  to  have  been  tilted  (into  the  xy  plane).  The  tilting 
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operation  is  specified  by  the  tensor: 

4 — ¥ 


0  0  ‘ 
cos  0  sin  0 
—  sin  0  cos  0_ 


(3.1) 


To  determine  the  effect  of  this  tilting  operation,  it  is  simplest  to  begin  with  the  resistivity 
tensor  equation  (2.21b).  Wc  inserted  a  subscript  p  to  equations  (2.21b),  (2.22)  and  (2.23) 
to  remind  readers  that  the  corresponding  matrices  have  been  specified  in  terms  of  the 
principal  coordinate  system  of  each  ellipsoid.  Let  us  now  compute  a  new  resistivity 
tensor  for  the  tilted  ellipsoid  as: 


fi(v)  _  <n  t{(v)  m7  _ 

Pr  —  -°0'P n  •  — 


S^e 


1  Mv)  —ZyH^/c  ZJlf/c 

ZvH[:'/c  i  /ip  -ZJl^/c+l/it 

L-zvh‘.7<-  zvh^/c  +  i/u 


,  (3.2a) 


where  the  subscript  t  refers  to  the  expression  that  is  found  subsequent  to  the  tilting  opera' 
tion.  In  the  above,  wc  have  introduced  the  definitions: 


1 /u!;1  s  l/»7 

1  /iiy  ]  =  (cos2  0)/ityJ  +  (sin2 
1  /uiv>  =  (cos2  0)/u{p  +  (sin2  0)1  iffi 

J/“  =  (^i  -  sil1 0cos0- 


H?  m 

H'P  =  up  cos  0  +  //(':>  Sin  0 
IW  =  -  sin  0  +  H'y  cos  0 

(3.2b) 


By  straightforward  inversion  of  equation  (3.2a),  one  obtains  the  conductivity  tensor 


associated  with  the  vth  tilted  ellipsoid  as 

’  1  1 

+("?Y 

ZJW 

zji';1  nyu? 
+  - 

z,h(:>  h^h':'  " 

rr 

'  c  / 

tt.c 

uc  cz 

»,c 

tic  c1 

«<')  -  — 

-ZXH[" 

ZXH';'  1 

[ 

z*n':} 

i 

,(3.3a) 

f  K* 

ti.C 

"1  2 
tiC  C 

uxu: 

+  \  c) 

uxc 

It  Xu  +  c2 

ZJi™  L  H'j'H';' 

1 1 

[(”  J 

1 

+  /HI")' 

uyc 

uc  c2 

tlxC 

ltx1l  c2 

nxuy 

\  V  ' 

where 

-r 

III 

> 

ytt:  —  (In 

i)(l>,i/  - 

(3.3b) 

and 

A(v)  a  1  +  uyu. 

:W?/c)2 

+  uxu:{H(yv)  c)2  -1-  iiAi/: 

(3.3c) 

According  to  the  model  adopted  here,  the  tilting  operation  is  unnecessary  for  ellipsoid 
v  =  4  in  the  valence  band;  its  appropriate  conductivity  tensor  is  given  by  equation  (2.16). 
In  principle,  the  various  mobilities  in  equation  (3.3)  should  be  indexed  with  a  superscript  v 
as  well,  but  the  latter  is  temporarily  omitted  for  the  sake  of  clarity. 

The  Scebeck  tensor  for  the  electrons,  equation  (2.22),  which  has  been  constructed 
subject  to  the  assumption  that  the  r?:  are  independent  of  r :,  is  isotropic  in  the  present 
approximation,  and  therefore  remains  unaffected  by  the  tilting  operation.  Also,  the  elec¬ 
tronic  contribution  to  the  thermal  conductivity  associated  with  the  charge  carriers  in  a 
particular  ellipsoid  can  be  obtained  via  equation  (2.23). 
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Before  summing  the  individual  into  a  total  conductivity  tensor,  it  is  necessary  to 
rotate  ellipsoids  v  =  2,3  so  that  their  principal  axes  coincide  with  those  of  the  crystal 
axes  (the  ellipsoid  v  =  1  is  already  properly  aligned).  For  this  purpose,  we  introduce  the 
transformations  (rotations  about  the  z  axis  by  120°  and  240°): 


(3.4) 


j(3)  ^(3) 

Their  effect  is  found  by  carrying  out  the  similarity  transformations  a  2  =  ®240  :  g(2  ®i2o 
on  equation  (3.3a).  Wc  then  obtain: 


1 

~  2 

Vi 

o' 

1 

—  T 

i 

K»|U» 

O 

11 

© 

r\ 

I® 

“Vi 

1 

2 

0 

®240  = 

Vi 

-k 

0 

.  o 

0 

1. 

-  0 

0 

1. 

120 


240 


‘V 

<1  2  -  -r- 


3  1  It 

1  +  A  + 

4  uxu.  4  nriis 

//:  1  I 

+  ~T  ~  A  > 
r"  4  ir 


i-v-1  1  i 

4  nxii:  4  wvu_ 

-Zv  If 


u: 

-/ 


-f-  — *—  -  ± 


s/3  J_ 

4  i< 2 


tfx  Uv\ 

2a  l+V-  c  “77 

2  ux  \  c  cl 

l  Zy  /  U  X  ll  A 

-A  (T  r - -)  + 

4  mv  \  c  C  J 

II XH,  ,31  1  Zv  //. 

+  *,--+■  —  -  — 
c2  2  2  ii  c 


+i2_ LtV3_L+ 

4  «xu.  4  uyu, 

Zv  H  HXIIV 
+——  +  -~± 


'.■4-up)- 

iAu.  4  nyii.  \  c  / 


4» 


3 

~4ir 


4«" 


/  -  //v  //v\ 

.(• 

J3ZV/  //x  //v\ 

±r~hv3—  -  ’  + 

4  i/v  \  c  c  / 

HJL  v  3  Zv  II.  1  1 

+-V  +  V— —  + . 

c  2  ii  c  2  hxh 


4  ii/  c  c  / 

1  ZV  Hx  H\ 

+i-b>/3— )+ 

4i//  c  c  / 

HXH.  ZVH .  J3  1 

+  — +  .  —  ±  ■ 


2w  c 


1  ZJ-HX  II A 

-T—  —  ±V3— It 

4  a,  \  c  c  / 


_  J3  Z,  H.  1  1 
2  uu, 


+-r"iT^--  + 

2  a  c 


1  +  (-f 

uxiiy  \  c  / 


in  which: 


1 


t//, 


1/ylY-  + 


+ 


uvir  uc  \  2  c  c ) 


(3.5b) 


In  the  above,  we  have  also  transformed  the  magnetic  held  components  which  were 
originally  expressed  in  terms  of  the  principal  coordinate  system  of  the  tilted  ellipsoids. 
The  magnetic  held  with  components  H,(a  =  .\\  y,  z)  is  now  specified  in  the  common 
coordinate  system  introduced  earlier.  In  carrying  out  the  transformation,  use  has  been 


(3.5a) 
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made  of  the  relationships 


one  should  also  note  that  H ^4)  =  H{P  =  H^k  =  .v ,y,z). 

In  principle,  we  are  now  in  a  position  to  obtain  the  total  a*  tensor  via  equation  (2.27). 
In  practice,  due  to  the  denominators  which  differ  for  various  v,  the  requisite  algebraic 
manipulations  beeome  excessive;  wc  must  resort  to  an  examination  of  speeial  cases  in 
order  to  obtain  tractable  results. 


4.  TRANSPORT  TENSORS  FOR  BISMUTH  WITH  H  —  kH 


In  this  seetion,  we  specialize  to  the  case  where  the  magnetic  field  is  aligned  with  the 
z  axis  of  the  crystal  coordinate  system.  With  this  specialization,  the  subsequent  mathe¬ 
matical  operations  are  considerably  reduced. 


We  first  demonstrate  that  for  the  conduction  band,  the  transport  integral  S(^  may 


be  replaced  with  ne/3,  where  nc  is  the  charge  carrier  density  in  the  conduction  band. 
This  identity  may  be  established  by  setting  //  =  1  in  equation  (2.13),  introducing  the 
changes  in  variable  =  z/kT,  i\B  =  and  integrating  by  parts;  one  thus  obtains: 


(4.1) 


The  central  expression  will  be  recognized  as  the  density  of  charge  carriers  in  a  valley, 
n(v\  Since  we  assume  that  there  are  three  equivalent  ellipsoids  in  the  conduction  band, 
it  follows  that  S(^  =  nJ3,  as  was  to  be  proved. 

In  now  adding  the  partial  conductivities,  we  return  to  equation  (3.3)  for  ellipsoid  1, 
to  equation  (3.5)  for  ellipsoids  2  and  3,  and  to  equation  (2.16)  for  ellipsoid  4.  In  these 
tensors,  we  set  II x  =  Hx  =  0.  Next,  we  restrict  ourselves  to  the  case  where  nc  =  nh  =  //, 


which  applies  to  pure  bismuth.  Because  of  the  considerable  overlap  between  the  valence 


and  conduction  bands,  this  equality  is  almost  always  experimentally  encountered.  Finally, 
we  utilize  the  fact  that  u{"]  g>  and  as  has  been  experimentally  well  docu¬ 

mented  [5,  16]. 

For  later  use,  we  now'  construct,  subject  to  the  restrictions  mentioned  above,  the  partial 
conductivities  of  carriers  in  the  two  bands.  These  quantities  arc  given  by  [see  also  equa¬ 
tion  (2.27)]: 


and  =^), 


v=  1 


The  requisite  algebraie  manipulations  are  straightforward,  though  rather  tedious,  and 
lead  to  the  following  results: 


-{Hz/c)a^ii(yn) 
Atf*/ 2 
0 


0 

0 


(4.2a) 


0 


!if[  1  +  ^hi^H;/c2l 


with 


D:  =  A  +  v^W/c2, 


(4.2b) 
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and 

As  1  -  u\n)u^/u(n)2. 

The  corresponding  tensor  for  holes  is  given  by: 

(Hjc)iW2  0  ' 

<4P)  0  , 

0  u[p)  A*. 

with 

A.  =  1  +  u(p)2H2Jc2. 

From  the  above,  one  can  then  construct  the  total  conductivity  tensor  as: 


a«\Hz)  = 


ne 


i,(p) 


(H,/c)u(p)1 

0 


(4.2c) 


(4.3a) 


(4.3b) 


'AA.u^ft  +  Dm <f> 
(Hjc)[ Ad?uf  - 


0 


-uijc)[Am<;%"'  -  d.Ap)1) 

AAm^/2  + 

0 


0 

0 


A,4n,[l  + 


uf'l %'Hh 

—?  J 


+  AtDtu[f) 


(4.4) 


The  resistivity  tensor  is  thus  given  by : 


Pz 


AAM{xn)/2  +  DmP 

-  Dm{p)2] 

0 


(HJc)IAm^  -  Dm'*} 
AAm^/2  +  Dm(p) 

0 


0 

0 

m2/D'x 


(4.5a) 


in  which : 

D’z  =  A,{wf[l  +  u^ufHl/c-]  +  D.u{p))D"_ ;  (4.5b) 

and 

D:  =  [AAX",/2  +  D:u[p)]2  +  [(H./c)(A;M*c",«‘r)  -  D.u(p)2)]2.  (4.5c) 

The  generalized  Seebeck  tensor  may  now  be  determined  according  to  the  relation 
pi  =  pi. :  (ole)  ■  Ve  +  P*)»  which  is  a  reformulation  of  equation  (2.28)  as  applied  to 

the  present  case.  The  partial  conductivity  tensors  are  given  by  equations  (4.2)  and  (4.3), 
the  total  resistivity  tensor,  by  equation  (4.4),  and  the  quantities  1pe  h  =  pe,/,|l]|  are  specified 
by  equation  (2.22).  On  carrying  out  the  indicated  operations,  one  obtains : 


where 


P;(//J  = 


d : 


PJH:) 

P„(-H_) 

0 


PXy(Hz) 

PJH:) 

0 


0  " 
0 

P-JH;l 


Vxx(H:)  =  [AVeA:u™/2  +  phD2u<»)[AAsu™/2  +  D.u^]+ 


+(HJc)2[peAzu*¥? 


-  phD:u(p)2] [ A^l4^)u'^,  -  DXP>2]; 


(4.6a) 


(4.6b) 


— 
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P =  (H;/c){[Ap A«?V2  +  -  Dt  itf]- 

+  Dzuip)][vA,t/;V;>  -  v„dAp)1]}; 

V„(HZ)  =  (D:7D:){p,A;i/W[l  +  ufufHl/c1]  +  p„D.  A^*}. 
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(4.6c) 

(4.6d) 


The  above  results  are  too  cumbersome  to  be  of  use  in  the  calculations.  Accordingly, 
we  consider  two  limiting  eases,  (a)  (u^Hjc)1  1  (“weak  field  ease’’),  and  (b)  (u^HJc)2  >  1 
(“strong  field  case'5)  where  /,  v  =  x,yfz  and  b  =  n,p.  One  should  remember  that  in  the 
latter  case,  all  quantum  oscillatory  phenomena  are  left  out  of  account,  so  that  the  strong 
field  results  should  be  regarded  as  approximations  to  the  region  of  intermediate  field 
strengths  where  quantum  effeets  are  not  yet  appreciable. 

A  limiting  proeess  in  which  terms  of  order  H 2  or  higher  are  neglected  leads  to  the 
results  shown  below : 


m 


i 


r  *  * 

mSo  *  ~  Ane(u(£)j2  +  u^)2 


~A{ufI 2  +  u[p))  (HJc)(u{?u{p  -  Au(p}2) 

-(II.Jc^uW  -  Au(f~)  A(t/;]f2  +  u Sf>) 


0 


0 


1 


•  «-  * 

«™o  Pl  “  2  +  i#’)2 " 

+  «Sf’)  -  (HJc)Axt;ht>\p,  -  ft)  x 


X  WJc)AxPi/jto'-  ft)x 
x  (i f?  +  AUt/l) 

0 


x  «>  +  Au^/2) 

A\ p.u?’/2  +  w&tben  +  <#>) 


a2(i/;]/2  +  t/'Y 


uj*’  +  Aif?  - 


;  (4.7) 


.42(ii,”)/2  +  +  p^/tuT’1) 

u'"'  +  /tut'” 


.  (4.8) 


Likewise,  on  retaining  only  the  highest  even  and  odd  powered  terms  in  (4.5)  and  (4.6), 
the  strong  field  results  read : 


lien  p.  =  — 
//,-x  '  ne 


2u'xf>ul;\HJc)1  4 uJftuf'nW  _  Aufl)(HJc) 

Au'?  +  2uf  ~  +  2u'*Y 

-4u'")«)<)  -  Au'f  )(UJc)  2uf;)uf(HJcf 


u(;\Au{xp,  +  2u(;)2 


u'”>  +  u<?>J 


(4.9) 


lim  p. 

//,-»0O 


r4<KX»  -  Ait»)(Ve  -  a) 

-2u(/,H?,(ft  -  ft)(W./c) 

0 

u(;\Au^<  +  2<>)2 

Au<?  +  2i#> 

2»Sf)«!r’(pe  -  ft)(H,A0 

4«<">(u?>u},>  -  M^Mft  -  Pd 

0 

(4.10) 

All'’'  +  2uf 

+  2UW)2 

0 

0 

vAn) 

+  w4rt 

< 

+  ui"  _ 

2 
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Attention  is  directed  to  the  faet  that  one  must  proeccd  via  the  exact  tensor,  equation 
(4.5),  before  taking  the  limit  Hz  go  beeause  in  several  instances  the  leading  term  of 
the  high  field  approximation  vanishes  identically.  In  these  situations,  the  corresponding 
entries  involve  terms  of  the  next  lower  order  in  Hz. 

A  diseussion  of  the  above  results  is  deferred  to  Part  II  of  this  paper. 


5.  TRANSPORT  TENSORS  FOR  BISMUTH  FOR  II  =  \HX  AND  II  =  j Hy 

Two  other  cases  which  remain  relatively  tractable  are  those  for  which  the  magnetic 
field  points  along  the  .v  or  y  directions  of  the  erystal  system.  The  procedure  is  exaetly 
the  same  as  outlined  in  Seetion  4  and  reviewed  again  in  Section  6,  cxeept  that  different 
magnetic  field  components  arc  set  equal  to  zero.  We  shall  therefore  only  eite  the  final 
results.  Despite  all  the  simplifications  used  so  far,  the  analogues  of  equations  (4.5)  and 
(46)  for  the  quantities  K\>  IV tire  extremely  cumbersome;  there  seems  little  point 
in  setting  them  down.  Instead,  we  proeccd  directly  to  the  weak  and  strong  field  limits 
eited  below.  We  find : 


lini  v.x  =  —7,0 

iix-o  ncD 


'2(11?'  +  Ai&)  0 

0  2(u‘">  +  Au'f) 

0  -(tifuf  -  2.4iifii<'’')//v. 


(tiftif  - 
..!(»«">  +  2, if) 


(5.1a) 


with 


and 


D"  =.  (tif  +  2nf)(nf  +  Aufy, 


(5.1b) 


lim  =  - - 

//*-*«  *  neDa 


3D* 

0 

0 


iifiif'Piif  +  tWWIJcf 
-(tfv/c)(3iifiif  -  At/r,tit,’>)+ 


(«x/c)(3hJ,)h?)  -  ^iif  nf )+ 
«<">M<f»[(2/t  +  1  )(A">  +  3iif](W,  (•)’ 


.  (5.2a) 


with 


D“  s  (3af  +  ufMiif  +  iif)  -  (ii^iifi'ify-ftfii'"’  +  iif); 


lim  (?,.  -  -  .. 

j/v-o  1  ‘  ncD 


2(iif  +  /tiif) 


-  tifuf] 


0 

2(iif  +  /tut'’') 
0 


-  (II y  <  )(2/4  wf  uf> -  (if  11?'] 
0 

/l(tif  +  2iif ) 


(5.2b) 


(5.3) 


lirn 

«x~*0 


lirn 


lim 
//*-.  o 


— 


lim  p„  = 

/fy-.so  *  n eD* 


Galvano-Thermomagnetic  Phenomena 

it<%(/>(3u,;>  +  u^iHyjc)2  -  {Hjcyu'? -  M- 

-(u?> +  3«<",Xi4">m?)  - 


-  3^u<f>u«)] 


-  (/yc)[u<")K<['’l«?)(  1  -  /l)  -  [2,-1  U^lWlW  - 

>*(|(?>  +  3»W) 


-  t««  +  3u<»w>„<»>  -  +  •*’>? 

'■  x  T  y  Ax  “*  („(«)  +  3H(m 


ul">  +  3ujrt 


This  disposes  of  (he  various  resistivities.  The  Sccbcck-Nernst  tensors  in  these  approxi¬ 
mations  are  given  by : 


o 


pA  - 


(ft"1 +  AuMp.u?  +  2pAu<f>)  o 

0  (ujr>  +  AuMp, *fi»  +  2pk«i'>)  -  +  2«</,Xp«  -  P h)HJc 

0  u?,uSf,(u'">  +  /tw^Kp.  -  p»)/4'c  (u?1  +  2«<f)XpX'"  +  ^P*u<rt) 


(5.5) 


P*  = 


[p«»?’  +  3p*u!/’,]0’t 


u?>  +  3  u<?' 
0 


0 


pei/«[/tu«  +  (2/1  +  1)1/'"’]  + 

+  PaI^Ou™  +  «<f>)+ 

+  {HJct  P«  -  p»)«5',l4",uS,,l4rt/«,"> 

—  (p«  —  p*)i/J!,,»Sc,’,{i((",/m,") 

-  (//,v/c)[3»l">  +  (2A  +  !)««]} 


V  Di 


(«?> + /t«irtxp.«?’  +  2p*«r) 


(«<">  +  /l4'’,XPe»?)  +  2PAU<f») 


-  +  /Ik^Xp,  -  \\)ujc 


0 

-  (p,  -  PA)«(",i/l'”{«'r)/u(">  + 

+  (HJcW?  +  3i(J")l 

pell,f(u,")  +  Alt/*)  +  PAWy^u'"’  + 
+  (2/1  +  1X4'>]  - 

-  (P.  -  Pa 

UW|((P){|(W  +  2uW)  x' 

X  (P,  -  Pa)^,/c 
0 

(«(;>  +  2»i?,xpti4”  + 

+  A  PauI1’1) 


(5.6) 


(5.7) 
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,  -  1 
lim  p  =  — 

ry  D 


{u(!Xu(J"  4  3wJ,))(3pci4p)  4  p*hJ») 


-  vhWJc) 


4  4P)[rX")(^Mi/,)  +  2/H")  +  4 


4-  pXgKjM  4  1)<  +}Au(xp)]][ 

u™  +  3  u(xp) 

-  2u™u{*\u™  4  Au(f')(pe  -  pj 
(H„/c)h(")(mJ|i)  4  3u<p)) 


PewSWO  -  /i)  4  mm?*  4  »im))  x 

x  (3pX")  4  p huT) 


2  »!cH)»y(Pe  ~  P.X//y/c)p/ 
(M?,4  3»Sf)) 


{-  -  p„X2Mp,»=/” 

(3uip)  4  i£V* 


ui"^(3wjn)  4  i/!p))  x 
x  (p<  -  VtMHy/c) 


~  pj 


{(3h<",-4  U(XP)XHP>  +  »xn))  X 
x  (p<fM?)  4  phAu(p>) 

4  <>u<p)(l  -  A)[2peu(f> 

4  pMn)  ±  u(xp)m 

(ii(?  4  3mJp))’ 


6.  PRELIMINARY  DISCUSSION 

The  above  derivations  provide  all  the  results  needed  in  the  derivation  of  appropriate 
figures  of  merit  for  Bi;  this  extension  of  the  work  and  the  comparison  between  theory 
and  experiment  is  deferred  to  Part  II  of  this  paper.  However,  at  this  stage  it  is  worthwhile 
to  review  in  outline  form  the  general  procedure  and  assumptions  used  in  the  derivations 
so  far. 

Starting  with  the  conventional  formulation  of  the  electric  and  thermal  currents  in 
terms  of  transport  integrals,  the  distribution  function  was  determined  by  solving  the 
Boltzmann  transport  equation  in  the  relaxation  time  formalism.  Allowance  was  made  for 
anisotropy  of  the  medium  and  for  the  joint  action  of  electric,  magnetic,  and  temperature 
fields.  The  results  so  obtained  specify  the  fluxes  J  and  in  terms  of  the  “forces”  V(J/e) 
and  VT ;  after  converting  from  to  Js  one  obtains  a  set  of  phenomenological  equations 
satisfying  the  Casimir-Onsager  reciprocity  conditions.  On  partial  inversion  of  these 
relations,  whereby  V(J/e)  and  Js  were  assigned  the  status  of  dependent  variables,  the 
phenomenological  equations  were  cast  in  a  form  that  permitted  the  identification  of 
the  resistivity,  Seebeck,  and  thermal  conductivity  tensors  in  terms  of  transport  integrals. 
Expressions  so  obtained  were  then  properly  summed  to  take  into  account  the  contribu¬ 
tions  arising  from  carriers  in  various  conduction  and  valence  band  ellipsoids,  allowance 
being  made  for  the  tilt  of  the  former  set  about  their  respective  binary  axes. 

The  final  results  are  based  on  the  following  simplifying  assumptions:  (a)  i:  is  assumed 
to  vary  quadratically  with  the  wave  number  vector  components,  (b)  only  contributions 
due  to  the  overlapping  valence  and  conduction  bands  were  considered,  (c)  the  relaxation 
time  formalism  was  used  in  solving  the  Boltzmann  equation,  (d)  all  phonon-charge 
carrier  interactions  were  presumed  to  be  accounted  for  in  specifying  the  relaxation  time, 
and  intervalley  scattering  effects  were  ignored,  (e)  hot  electron  effects  and  quantum 
effects  in  the  high  magnetic  field  region  were  not  included.  To  obtain  tractable  analytic 
expressions,  we  also  (f)  considered  only  the  special  case  in  which  the  magnetic  field  is 
aligned  with  each  of  the  crystal  symmetry  axes,  (g)  neglected,  where  permissible,  1/hJ* 
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relative  to  l/w*w)  and  set  u ^  =  vff\  where  the  uib)  are  electron  or  hole  mobilities  for 
H  =  0  along  the  indicated  crystal  symmetry  axes,  (h)  restricted  ourselves  to  the  intrinsic 
case  by  setting  n  =  p ,  and  (i)  specialized  to  the  case  of  low  and  high  magnetic  fields. 
Finally,  (j)  for  the  weak  field  limiting  case,  only  terms  of  order  H°  and  Hl  were  retained. 
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Resum 6— On  a  appliqu£  au  bismuth  le  formalisme  du  temps  de  relaxation;  les  r&ultats  sont  utilises,  dans 
la  partie  II,  pour  determiner  le  facteur  de  m£rite  du  bismuth  dans  les  processus  de  conversion  d’^nergic. 
En  utilisant,  pour  le  bismuth,  le  module  de  Jones-Shoenberg,  on  derive  des  expressions  analytiques  de  la 
resistivity  yiectrique,  de  la  conductivity  thermique,  et  des  coefficients  de  Hall,  de  Seebeck  et  de  Nemst. 
On  rysout  liquation  de  transport  de  Boltzmann  donnant  la  fonction  de  distribution  perturbee,  en  utilisant 
des  temps  de  relaxation  anisotropes.  Le  rysultat  est  alors  introduit  dans  les  intygrales  de  transport  repry- 
sentant  le  courant  yiectrique  et  le  flux  d’ynergie,  pour  aboutir  aux  equations  phenomenologiques  valables 
pour  chaque  groupe  de  porteurs  de  charge  associe  k  un  ellipsolde  donny.  On  additionne  les  contributions  de 
chaque  groupe  de  porteurs,  dans  le  systeme  de  coordonnees  du  cristal,  pour  obtenir  les  coefficients  de  trans¬ 
port  mentionnys  ci-dessus.  Pour  obtenir  des  expressions  analytiques,  il  est  necessaire  de  considerer  Ics  cas 
particulars  oCi  le  champ  magnetique  est  aligne  avec  chacun  des  trois  axes  de  symetrie,  et  de  passer  k  la 
limite  des  champs  magnytiques  trys  faibles  ou  trys  forts. 

Zusammenfassung  Die  Transport-Theorie,  begriindet  auf  dem  Relaxationszeit-Formalismus  wird  auf 
Wismut  angewendet.  Diese  Erbebnisse  werden  dann  in  Teil  II  benutz,  um  die  Giiteziffer  von  Bi  bei  Energi- 
cumwandlungs-Prozessen  zu  bestimmen.  Unter  Benutzung  des  Jones-Shoenberg-Models  fiir  Wismut 
werden  analytische  Ausdrlicke  fur  den  spezifischen  elektrischen  Widerstand,  die  Warmeleitfahigkeit  und 
fur  die  Hall-,  Seebeck-  und  Nernst-Koeffizienten  abgeleitet.  Die  Boltzmann’sche  Transportgleichung  wird 
gelost  fiir  gestorte  Verteilungsfunktionen  unter  Benutzung  von  anisotropen  Relaxationszeiten.  Dieses 
Resultat  wird  dann  in  die  Transport-Integrale  fur  den  elektrischen  Strom  und  den  Energiefluss  eingefuhrt, 
um  phanomenologische  Gleichungen  fur  jede  Art  von  Ladungstragem  zu  erhalten,  die  mit  einem  gegebenen 
Ellipsoid  verbunden  sind.  Die  Beitrage  jeder  Gruppe  von  Ladungstragem  werden  dann  addiert  in  dem 
gemeinsamen  Symmetrie-Koordinatensystem  des  Kristalls,  um  die  oben  erwahnten  Transportkoeffizienten 
zu  erhalten.  Um  analytische  Ausdriicke  zu  erhalten,  war  es  notwendig,  die  Sonderfalle  durchzuarbeiten, 
in  denen  das  Magnetfeld  mit  jeder  der  drei  Symmetrie-Achsen  ubereinstimmt  und  tiberzugehen  auf  die 
Grenzfalle  sehr  niedriger  oder  sehr  hoher  Magnetfelder. 


